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Abstract 

Proper tuning of the orbital characteristics of the three spacecrafts that constitute the usual 
triangular configuration of the space-borne gravitational-wave detector LISA, could minimize the 
breathing mode of its arm-lengths. Since the three spacecrafts form three pairs of interferometric 
arms, we have the freedom to minimize whichever combination of arm-length variations that might 
be useful in signal analysis. Thus for any kind of time delay interferometry (TDI), that is chosen 
to be used in analysing the data, the optimal orbital characteristics could be chosen accordingly, 
so as to enhance the performance of the gravitational wave detector. 

PACS numbers: 04.80.Nn 
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I. INTRODUCTION 



The joint ESA-NASA future mission to launch a spaceborne gravitational wave antenna, 
known as LISA (Laser Interferometer Space Antenna), is expected to offer us invaluable in- 
formation for the Universe. The detector will be able to monitor low frequency gravitational 
waves, the source of which could be either supermassive black hole binaries at cosmological 
distances [ll, Q], white dwarf and/or neutron star binaries [3], or primordial waves of cosmo- 
logical origin y,[5|. Detection of such gravitational waves could be achieved by placing three 
spacecrafts into three distinct Earth-like orbits, so that they form an equilateral triangle of 
almost constant size Q|, and monitoring interferometrically the tiny variations of distances 
between any pair of such spacecrafts that are induced by gravitational waves passing by 
LISA j?]. Since the frequency band that this detector is sensitive at is in the region 10~^ to 
10^^ Hz, any possible relative motions of spacecrafts that occur at much lower frequencies 
of order ~ 10"'' Hz (due to orbital periodicities) could be filtered out from the signal [s], 1^. 
However, despite the fact that the large beat modes due to Doppler shifts caused by arm- 
length variations could be effectively eliminated, the same arm-length variations could cause 
unpleasant complications while attempting to reduce internal noise by various schemes of 
time delay interferometry (known as TDI's) [ll| . This happens because, while the laser 
beams travel back and forth across the arms, they pass through each spacecraft at different 
times due to arm-length changes; thus any reference frequency variations of local lasers do 
not exactly cancel out in the corresponding combination of signals, as it would happen if 
there was no flexing of the arms. Consequently, if the variability of the size of the triangular 
formation could be minimized by suitable adjustment of their orbital characteristics, this 
would lead to a suppression of the internal noise of LISA, which, in its turn, would be highly 
beneficial for its performance as a gravitational wave detector. 

We start by writing down the position of each spacecraft, accurate up to second order 
with respect to its eccentricity. Each spacecraft moves on a slightly elliptical orbit with 
semi-major axis equal to 1 AU, the plane of which is slightly inclined with respect to the 
plane of the ecliptic. Finally the orbit of each spacecraft is rotated by an angle of ±27r/3, 
with respect to the other two orbits, on the ecliptic plane. The initial position of the 
three spacecrafts on their corresponding orbits is such that they form an equilateral triangle 
that remains equilateral to first order with respect to the eccentricity e of the orbits. The 
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proposed configuration of LISA, tliat is widely used in corresponding analyses [g], assumes 
that the plane of the triangular configuration forms an exact 60° angle with respect to the 
ecliptic plane, with the triangle being exactly equilateral initially. Indeed this configuration 
ensures stable distances between spacecrafts to first order with respect to the eccentricity 
of the orbits e. However, this configuration leads to a breathing mode of the arms with 
amplitude of order 12|. This arm flexing end up generating noise, by one way or another. 



in detector's output 



It is to verify that each spacecraft has 3 extra degrees of freedom that could be used 
to minimize whichever arm-length variation might one choose. Of course this fine tuning 
of the orbital characteristics should be one order of e higher than the initially proposed 
value of the characteristics themselves, so that the invariance of the arm-lengths to order 
e is not destroyed. These nine, altogether, degrees of freedoms could be chosen to be the 
eccentricities of the orbits, their inclinations, and the initial angle position of each spacecraft 
along its own slightly eccentric and inclined orbit. The choice could be such that the initial 
configuration deviates from being equilateral to order e, and/or its plane inclination deviates 
from the 60° to order e, as well. On the other hand such an adjustment could reduce the 
time-variation of any arm-length, or the relative variation of any pair of arm-lengths, at its 
minimum value. 

The rest of the paper is organized as follows. In Sec. [IT] we repeat the basic calculations, 
in a form that could be later extended to higher order with respect to e, that show why the 
initially designed configuration leads to a time- independent equilateral triangular configu- 
ration to order e. By expanding the spacecraft distances to order e^, we compute the time 
variation of the arm-lengths at this order. In Sec. IIIII we introduce nine extra parameters 
of order that modify the initial positions of the three spacecrafts. Then we calculate 
once again the distances between each pair of spacecrafts as a function of time, which are 
now parametrized by the six (out of nine) essential parameters that determine the initial 
location of the spacecrafts. Finally, in Sec. IIVI we show how we could optimally choose the 
fine tuning parameters so as to get the minimum contribution to noise from the variations 
of the arm-lengths, depending on the TDI scheme that one might choose to use in signal 
extraction. 
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II. THE INITIALLY DESIGNED ORBITS 

The orbit of each spacecraft is a Keplerian elhpse, if one ignores the gravitational at- 
traction due to planets (mainly Earth, and Jupiter) and relativistic effects. Assuming that 
all spacecrafts are moving on orbits that have semi-major axis of 1 AU, so that all return 
at their initial positions after one year, and are not drifting away secularly with respect to 
Earth, their distance from the Sun could be written as 



where a, the semi-major axis, is common for all spacecrafts and is equal to 1 AU. The 
eccentricity, e, is common for all spacecrafts, as well, for the proposed configuration and is 
the only parameter that uniquely characterizes all three orbits (see below). Finally, ^ is the 
so called eccentric anomaly, which is an angle parameter that determines the position of the 
spacecraft along its elliptical orbit. Now the orbital plane of each spacecraft is inclined with 
respect to the ecliptic plane by an angle Ae, where A is a number that will be determined 
by demanding the configuration of the three spacecrafts to be an equilateral triangle with 
invariant size to first order with respect to e. The three inclined orbital planes are rotated 
by an angle of ±27r/3, with respect to each other, on the ecliptic plane (see Fig. [T]). By using 
the eccentric anomaly to describe the position of each spacecraft along its orbit, instead of 
the polar angle 9, the calculations that lead to the distance between two spacecrafts as a 
function of time are made easier since by conservation of angular momentum the angular 
position ^ on a keplerian orbit as a function of time is given by 




(1) 





while by using the eccentric anomaly parameter ^ which is related to 6 via 



^ 1 — e cos 9 ' 



(3) 



the cartesian coordinates of the orbital position are 



x= r cos 6* = a(cos^ + e). 



(4) 



y = r sin 9 = aV 1 — sin ^, 



and the integral of Eq. ([2]) is easily computed to yield 
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(c.f. [14!]). The last expression is easier to use, than Eq. ([2]), to compute t(^), in contrast to 
t{6). The above transcendental equation could not be exactly inverted in a closed analytical 
form. However, this could be achieved in the form of a power expansion with respect to e, 
yielding 



^ = Lut + ^0 + {sm{iut + ,^0) - sin ^0) X 

(-e + cos{ujt + ^0) + 0(e^)) 



(6) 



The cartesian coordinates of each spacecraft on the heliocentric system, with the z-axis 
perpendicular to the ecliptic plane, are thus given by 



(7) 



where xip is the position of the i-th spacecraft on its orbital plane, that is 



/^«cos0w\ f a(e + cos^«) ^ 



op 



sin 9^^^ 




ay/1 — sin ^ W 




(8) 



where the expressions in the last array come from Eqs. (jl]). The matrix produces an 
inclination of the orbital plane by an angle Ae with respect to ecliptic plane around y-axis. 
Therefore, 



^ cos(Ae) sin(Ae) ^ 



(9) 



1 
^ — sin(Ae) cos(Ae) j 

Finally, i?rot is the matrix that rotates the orbit of the i-th spacecraft by $1 = 0, $2 = 27r/3, 
and $3 = 47r/3, respectively, on the ecliptic plane, that is 



R 



rot 



V 



(10) 



cos $j sin $j 
— sin $j cos $j 

V 

By combining the coordinates of each spacecraft given above, and choosing the initial angular 
positions ^q*'* (the superscript refers to the i-th spacecraft) to be 



ISO )S0 )S0 i 



271 VS in VS ^ 



(11) 



along with the common inchnation parameter A = \/3, we obtain the following time de- 
pending distance of each pair of spacecrafts to second order with respect to e: 




ruit) = 2^306 - (l9cos(cjt) - 19\/3sin(cjt) - 2cos(3cut)) ae^ + 0(e^), (12) 



r3i(t) = 2V3ae- — (l9cos{ujt) + 19V3sm{ujt) -2cos{3ujt))ae'^ + 0{e^). (14) 



The choice of initial parameters for the three spacecrafts, mentioned above, is the one that 
leads to the essential advantage of the configuration; that is to keep the distances equal, 
and time-invariant to lowest order (first order) with respect to parameter e. However, the 
arm-lengths vary with time to order e^. The magnitude of the arm- length oscillation is thus 
of order 10^ km as shown in Figure [2l 

III. LOWERING THE AMPLITUDE OF ARM FLEXING 

In this section we will use our freedom of choosing the initial positioning of the three 
spacecrafts, in order to minimize whichever distance variation we might like. Of course by 
placing two spacecrafts on a circular orbit around the Sun the distance between them will 
remain fixed, but in order to sense the quadrupole nature of a gravitational wave we need 
al least one more spacecraft that is not along the same line of the former two ones. Thus, 
by placing a third spacecraft in an orbit that is inclined with respect to the orbit the other 
pair, we will have a time varying arm-length between the third and each one of the other two 
spacecrafts, the overall variation of which will be of the order of the initial arm-length. The 
clever symmetric configuration of the three slightly inclined and slightly non-circular orbits 
that was discussed in the previous section manages to keep all the arm- lengths constant, 
at least at the order of magnitude of the arm-lengths themselves. On the other hand by 
trying to achieve a very symmetric configuration we have ignored any possible freedom we 
still have to shift our orbits so as to reduce the variation of distances to even higher order. 

Speaking of freedom of initial positioning, the three elliptical orbits could be a little 
different, with respect to eccentricity and with respect to inclination, with each other. Of 
course the fine tuning of the corresponding six parameters will be of order e^, so as to keep 



r23(t) = 2V3ae + — {7cos{ujt) + cos{3ujt))ae'^ + 0{e^), 



(13) 
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FIG. 1: The figure depicts the orbits of the three spacecrafts labeled 1, 2, and 3 respectively. The 
triangular configuration is shown at an instance when spacecraft 1 is at its aphelion (^i = 0) while 
the other spacecrafts are located at angles 6*2 = 27r/3 + ©(e), and 63 = Att/3 + 0{e), respectively. 
The three elliptical orbits are inclined with respect to the ecliptic plane (shown in the figure); thus 
the part of the ellipses that lie above the ecliptic plane are drawn thicker than the parts that lie 
below. Also the three orbits are rotated with respect to each other by exactly ib27r/3, as shown in 
the figure. The inclinations and the eccentricities presented in this figure are highly exaggerated 
to clearly illustrate the geometry of the orbits. 



the main characteristic of the configuration; namely the constant value of the arm-lengths 
at order e. Furthermore one could also loosen the exact symmetric placement of the orbits 
by a rotating angle of ±27r/3 with respect to each other. However these angles could be 
kept invariant, and we could alternatively adjust the initial angular positions ^q*^ to order 
e\ 

Henceforth we will assume that the three orbits are characterized by the following orbital 
parameters: 

ei = e, e2 = e + Q;2e^ 63 = e + age^, 

(15) 

\i = V3 + Ae, A2 = 73 + (326, A3 = 73 + P^e. 
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FIG. 2: This diagram shows the variation of the three arm-lengths for the configuration of LISA 
that was initially designed. The range of arm-lengths' variation is 6.6 x 10^ km, for ri3, and ri2, and 
2.1 X 10^ km, for r23, respectively. The first two pairs have similar time evolution due to symmetric 
arrangements of the corresponding sides of the equilateral triangle. Actually ri2(t) = ri3(— t), if 
t = corresponds to configuration shown in Figure [2j 



We have actually chosen the eccentricity of the first spacecraft's orbit as a reference for 
the other two. This eccentricity will play the role of the small expansion parameter in all 
positional expressions. Next we define the initial positions of the spacecrafts along these 
slightly deformed, and differently inclined orbits by the initial angles 

= 0, 

ff = |-fe + o.e^ (16) 

where once again by setting ^q^"* equal to zero, we have just adjusted the initial time of the 
configuration to coincide with the aphelion of spacecraft 1. As we see now there are only 7 
parameters that we could tune to make the configuration comply with our demands. The 
assumed shifts of the orbital parameters of the three spacecrafts are independent to each 
other, and thus they exhaust the whole freedom we have to shift the initial positions of the 
spacecrafts. 

Actually the proposed configuration that forms initially an exact equilateral triangle in- 
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clined by 60° with respect to the equatorial orbit corresponds to some specific relations 
between these 7 new parameters since, as is shown in Figure [2], omission of all these param- 
eters does not lead to three equal arm-lengths at t = 0. 

Repeating once again the computations of the previous section with all seven new param- 
eters introduced in the formulae for Xop, R[^c^ ^o*'*' obtain the following expressions 
for the three distances as functions of time: 



( ri2{t) ^ 



2V3i 



ae 



1 



+ 



ae 
32 



A(t) + 0(6^ 



(17) 



with 



A{t) = Co + Ci cos(cjt) + C2 cos(2cjt) + C3 cos(3cjt) + Si sin(cut) + S2 sin(2o^t). (18) 

The 3x1 column matrices C^'s and S^'s, through which we have decomposed the time- 
depending second-order column matrix A(t), correspond to the following expressions that 
depend only on the seven fine-tuning parameters introduced above: 



Co = 12 



8a2 + 73 (/3i + (32) 
Sag + V3{f3i + Ps) 
y8{a2 + a3) + V^{p2 + P3) J 

C2 = I2V3 



Ci 

Pi 



-57-8(3a2 + 2v^72) 

-57 -8(3«3- 2^373) 
142 - 16(3(^2 + as) + 2^3(72 - 73)) I 



C, = 6 



Si 



\-il32+l33) J 

( 19y3 + 8V3«2 + I672 \ 
19^3 - sVSaa + I673 
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1 

V 

-Pi - 2Ps 
P2-P3 



(19) 



\ 



An obvious optimization choice of parameters is (3i = (32 = ((^3 = 0, since then the com- 
ponent of the variation of arms that oscillates at frequency 2u! vanishes (C2 = S2 = 0) 
without affecting the rest time- depending components Ci,C3,Si. On the other hand the 
only component of arm-length oscillation with frequency Su cannot be adjusted through 
suitable choice of the parameters. Finally the magnitude of the components corresponding 
to frequency uj could be adjusted by varying the value of specific combinations of 02, 72, and 
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as, 73- More specifically by using the following replacements 

X2 = 57 + 8(3a2 + 2y372), 



(20) 



;^3 = 57 + 8(3a3-2v^73), 

where the two new parameters X2,X3 could be adjusted independently to each other, the 
lowest-order arm-length variations oscillating at frequency u turn out to be 



ae 



\ 



ae 



32^3 



COs(co't)+ 



V " / 



sin(co't) 



(21) 



(22) 



-X2 

156 - (x2 + X3) J 

2x2 cos{ujt — 27r/3) 
2x3 cos(ci;t — 47r/3) 
\^ [156 - (x2 + Xs)] cos{ujt) J 

Therefore, all possible optimizations could be done by suitable adjustments of X2,X3 alone, 
that is, by simultaneous adjustments of a^'s and 7j's. For example, by choosing X2 = Xs = 0, 
two out of three arms are oscillating with the lowest possible amplitude, which is a mere 5% 
of the amplitude of the initially designed configuration. The price though is that then the 
third arm is oscillating with an amplitude that is ~ 3.4 times higher than the corresponding 
amplitude of the initial configuration. On the other hand if we need to use all arms in our 
signal analysis, we could make a compromise by a suitable choice of Xi's and manage to 
decrease the oscillation amplitude of all arms at an optimized ratio. In Figure [3], we have 
plotted the time- varying arm- lengths in a period of one year, for two choices of Xi's- 

IV. REDUCING THE NOISE IN VARIOUS TDI SCHEMES 



The optimal geometric configuration, that we could achieve by suitable initial positioning 
of the three spacecrafts, depends on the specific TDI (time delay interferometry) scheme that 
we decide to use in order to extract the signal from the internal noise of local lasers in each 
spacecraft. For example, if we had chosen to synthesize the phase-differences in an equal-arm 
four-link Michelson scheme interferometry (see 
mode of the arms would be 
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ll|) the noise induced by the breathing 



2C,{t){L2{t)-L3{t)), 
10 



(23) 




FIG. 3: This diagram shows the time varying length of the three arms when the orbital charac- 
teristics of the three spacecrafts are a little bit different than the ones of the initial design. The 
left diagram shows the extreme case X2 = X3 = Oi where two of the arms are oscillating at a small 
fraction of the oscillation amplitude of the initial configuration (only by 0.2 x 10^ km in amplitude), 
while the third arm is oscillating with an amplitude which is more than 3 times larger than the 
amplitude in the initial configuration. The right diagram shows a more balanced fine tuning of the 
parameters {x2 = Xs = 39) where all arms are oscillating with an amplitude ~ 60% lower than the 
maximum amplitudes of initial configuration. 



assuming that the spacecraft 1 is the corner station of the corresponding Michelson inter- 
ferometer. Li{t) is the time depending length of the i-th arm, while Ci{t) describes the 
time-dependent fluctuations of the i-th laser frequency (assuming there is only one reference 
laser in each spacecraft). The lasers of LISA mission are designed to have single-sided spec- 
tral density of order 30 Hz/v^Hz. The initial configuration has a time varying arm- length 
difference which has an amplitude of 3.3 x 10^ km (c.f. Figure [2]). On the other hand in the 
optimized configuration with X2 = Xs = although these two arms actually breathe, they 
have continuously exactly the same length (at least to order which corresponds to our 
approximations). The same holds good also for an unequal-arm eight-link Michelson com- 
bination, denoted X by Armstrong, Estabrook, and Tinto 15|, which is a second-generation 
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TDI scheme [Ul], since the antisymmetric combination L2L3 — L3L2 is again continuously 
zero (to the sanae order of approximation), and thus the internal laser noises cancel out 
(c.f. Eq. (12) of [11]). 

Let us examine one more second-generation TDI scheme that is, now, not symmetric with 
resp ect to a specific pair of arms. For example, the eight-link Relay scheme, denoted U in 
Iq . leads to the following TDI noise due to laser internal noise: 

C^it)[{Lr + U){Ly + L2O - {Ly + L2'){Lv + ^1) + ULv - U'Li]. (24) 

By expressing the three arm-lengths as Lj = Lj/ = Lq + li{t) where Lq is of order ae while 
li{t) is of order ae^ and the numbering is such that it corresponds to the i-th component of 
the column matrix of Eq. ( !22l) the above expression for the noise yields 



g3(t) ^^°^"' [(156 + X2- 2x3) sin^t + v^(2x2 + Xs) cosa;t] + O(e^). (25) 

To compute the above expression, only the Fourier components corresponding to frequency 
uj has been written down, since the Fourier component of arm breathing with frequency 3u; 
is the same for all arms (c.f. C3 of Eq. fll9p ). and the corresponding terms cancel out. From 
expression (l25l) it is easy to verify that by choosing X2 = —156/5 and xs = 312/5 we could 
nullify the noise of this TDI scheme to this order; namely a^e^. It should be noted that the 
specific choice of numbering of the arms with respect to the assumed distance Vij that we 
have used in our analysis leads to the specific optimizing parameter values X2,3 that we have 
found. 

Actually, all second generation eight-link TDI schemes lead to similar expressions for 
the noise, which could be written as antisymmetric products of arm-length variations and 
arm-lengths [ll]. These products could be expressed as a combination of sin cut and cos cut 
terms with corresponding factors that depend on the two parameters X2,3- Therefore there 
is always a suitable combination of the X2,3 parameters that eliminates the noise to that 
order, which means that by suitable fine initial positioning of the three spacecrafts we could 
depress the laser noise at the level of Ca^e'^uj. This is the best optimization we could achieve 
for a specific second generation eight-link TDI scheme based on the kinematics of LISA. 
As an order of magnitude, this means a reduction in the noise of LISA by e ^ 1/100 with 
respect to a non-optimizing positioning of the spacecrafts. 
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V. CONCLUSIONS 



In this short paper we have shown that we could adjust the orbital characteristics of the 
three spacecrafts which consist LISA detector at one order, with respect to e, higher than 
what is initially designed, in order to achieve specific kinematical properties. Namely, we 
could make the breathing mode of detector arms be optimized with respect to noise induced 
in the signal through any TDI scheme used to reduce the noise implications. We have shown, 
by presenting a few examples, that suitable initial positioning of the three spacecrafts could 
reduce the noise, due to lasers, by two orders of magnitude with respect to initial design. We 
should note though that since the positioning of spacecrafts could not be changed throughout 
mission's lifetime, only a specific TDI could be highly optimized. If another TDI scheme is 
used simultaneously to analyze some signal the benefits of the fine-tuned kinematics will not 
be equally highlighted. Hence the choice of kinematics should be based on the TDI scheme 
that will be most often used in signal analysis. 
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